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Abstract 

Ab initio electronic structure calculations are performed to support 
and to help interpret the experimental work reported in the proceeding 
manuscript. The CCSD(T) approach, in conjunction with a large ba- 
sis set, is used to compute spectroscopic constants for the X x S + and 
3 n states of ZnO and the A 2 S+ state of ZnO - . The spectroscopic 
constants, including the electron affinity, are in good agreement with 
experiment. The ZnO EA is significantly larger than that of O, thus 
relative to the atomic ground state asymptotes, ZnO has a larger Do 
than the 1 E + state, despite the fact that the extra electron goes into 
an antibonding orbital. The changes in spectroscopic constants can be 
understood in terms of the X 1 E + formally dissociating to Zn 'S + 0 
1 D , while the 3 IT and 2 E + states dissociate to Zn X S + O 3 P and Zn 
1 5 and O - 2 P, respectively. 

I. INTRODUCTION 

Fancher, de Clercq, Thomas, and Bowen 1 recently measured the photoelectron 
detachment spectra of ZnO - . Their electron affinity (EA) of ZnO is 2.088±0.010 eV, 
which is significantly larger than that of O (1.461 eV 2 ). This means that the D 0 
of ZnO - is 0.63 eV larger than that of ZnO (1.61T0.04 eV 3 ), if both are computed 
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with respect to the atomic ground states. While the D 0 was larger, the r e of ZnO - 
is 0.07 A longer than that of the neutral. Consistent with the longer bond length, 
their ZnO - u> e value (625 cm -1 ) is smaller than that for the neutral 4 (811 cm 4 ). The 
spectra also exhibited some features that might indicate a very low-lying electronic 
state of ZnO. 

Bauschlicher and Langhoff previously 3 studied ZnO and found the ground state 
to be 1 E + , however, the 3 II state was extremely low-lying. The Zn ls-Sd and 0 Is 
and 2s orbitals are not involved in the bonding, thus the valence occupation of these 
two states is 9cr 2 47r 4 and Q^lOcr'iw 3 , where the 9cr orbital is a bond between Zn 
4s and 0 2 per orbitals, the 10<r orbital is the antibonding combination of the same 
orbitals, and the 47 t orbital is mostly the 0 2 j>7t orbital with some donation to the 
Zn 4p7r orbital. Given these occupations for the two lowest states of ZnO, the ground 
state of ZnO - is 2 E + with an occupation of . Thus relative to the ground 

1 £+ state, ZnO" has the additional electron in the 10<r antibonding orbital. Given 
the occupations of the ground states of ZnO and ZnO - it appears strange that the 
D 0 of ZnO - should be significantly larger than that of ZnO. While the relative D 0 
values seem inconsistent with the occupations, the relative r e and w e values show the 
expected trends. 

In fight of the unexpected results, we have performed accurate calculations for 
ZnO and ZnO - to confirm the experimental observations and to understand the 
bonding. We also perform accurate calculations on both the 1 E + and 3 II states of 
ZnO to position these states and therefore help interpret the spectra. We also report 
results obtained using density functional theory (DFT), to compare the results of 
these inexpensive calculations with those obtained at the higher level of theory. 
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II. METHODS 


The DFT calculation are performed using the hybrid 6 B3LYP functional 7 in con- 
junction with the 6-311+G(2d) basis set 8-11 . These calculations are performed using 
Gaussian 94 12 . The accurate calculations use the coupled cluster singles and dou- 
bles approach 13 including the effect of unlinked triples, which is determined using 
perturbation theory 14 ; this approach is denoted CCSD(T). The restricted open-shell 
CCSD(T) approach 15 ’ 16 is used. In these calculations the O ls-like orbital and Zn 1 s- 
like to 3p-like orbitals are not correlated. In the CCSD(T) calculations, the Zn basis 
sets is the (20s 15p 9 d 6/ 4g)/[7s 6p 4 d 3/ 2 g] atomic natural orbital 17 (ANO) set 
that has been described previously 18 . The O basis set is the augmented correlation- 
consistent polarized- valence quadruple zeta (aug-cc-pVQZ) set developed by Dunning 
and co-workers 19,20 . Only the spherical harmonic components of the basis sets are 
used. These CCSD(T) calculations are performed using MOLPRO 96 21 . 

In order to gain insight into the bonding, the modified coupled pair functional 22 
(MCPF) natural orbitals are plotted. These calculations are performed at the B3LYP 
geometries. The basis set is that used in the B3LYP calculations with the diffuse Zn 
functions deleted and the two O 3 d functions replaced by one with an exponent of 0.8. 
Using this basis set, self-consistent-field (SCF) based constrained space orbital varia- 
tion 23 (CSOV) calculations are also performed. These MCPF and CSOV calculations 
were performed using the Molecule- Sweden program system 24 . 


III. RESULTS AND DISCUSSION 

The results of our calculations are summarized in Table I along with the previous 
CPF results 5 . The present calculations use a much larger basis set and the CCSD(T) 
approach yields a much better description of the electron correlation than does the 
CPF, thus the current CCSD(T) results are much superior to the old CPF results. 
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For the 1 E + state, which is not as well described by a single configuration as the 
3 II state, CCSD(T) yields a much larger D 0 , a higher frequency, and a shorter bond 
length than the CPF. The computed D 0 value is in excellent agreement with the 
experimental value 3 of 1.61T0.04 eV. We expect our computed value to be too small by 
approximately 0.1 eV, therefore, we conclude that the experimental value is probably 
slightly too small. Unlike Do, the computed u> e value (727 cm 1 ) is significantly 
smaller than the experimental value 4 of 811 cm -1 ; a difference of this magnitude is 
unexpected. The CCSD(T) value is significantly better than the older CPF value. The 
3 n state, which is better described by a single reference, shows smaller changes with 
the improvements in the calculations, and as a result, the T e value at the CCSD(T) 
level is much larger than reported previously. Thus while the old calculations suggest 
electron detachment should form two states at about the same energy, the current 
CCSD(T) results shows that the 3 II state is about 0.26 eV above the 1 S + state. 

The computed change in r e between ZnO - and ZnO (0.045 A) is smaller than 
found in experiment (0.07 A). In addition, the computed u> e of the 2 E + state (664 cm 
is only in reasonable agreement with experiment (625 cm 4 ). Unlike the Y 1 E + state, 
the computed value is too large for the 2 E + state. The norm of the singles amplitudes 
suggests that the 2 E + state is the most difficult to describe, and therefore probably 
has a slightly larger error than the two neutral states. The computed EA is too small, 
as expected. However part of this error arises from the error in the EA of the O atom 
(1.403 vs 1.461 eV). If we shift the Zn+O and Zn+O - asymptotes to agree with 
experiment, the corrected ZnO EA (2.03 eV) is in better agreement with experiment 
(2.088 eV). Thus the CCSD(T) calculations confirm the experimental observation 
that the EA of ZnO is significantly larger than that of O, that the ZnO - r e value is 
larger and the uj c value is smaller than those found for the ground state of ZnO. 

In Table I we also give results obtained at the B3LYP level. The r e values are in 
reasonable agreement with the CCSD(T), in fact for the 1 E + state the B3LYP results 
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agree better with the CCSD(T) than do the old CPF results. The accuracy of the 
uj c values is mixed, for the 1 E + state the agreement is very good, while for 3 II and 
2 E+ states the B3LYP results are about 70 cm” 1 smaller than the CCSD(T). The D 0 
values are too small and the B3LYP incorrectly places the 3 II state below the 1 E + 
state. The ZnO EA at the B3LYP level is too large, as is the EA of O atom. However, 
the corrected ZnO EA value agrees with experiment about as well as the CCSD(T) 
result. Overall the B3LYP results are reasonably accurate. 

The MCPF valence natural orbitals for the 3 n state of ZnO are plotted in Figs. 1, 
2 and 3. The orbitals for the a E + state of ZnO and for the 2 E + state of ZnO - are 
shown in Figs. 4-8. The MCPF populations are given in Table II. From the plots of 
the a orbitals it is clear that the Zn undergoes 4s4p hybridization, and these form 
a bonding and antibonding orbital with the O 2 per orbital. The 47T is mostly O 2jm 
with some donation to the Zn. The orbitals of ZnO - are somewhat more diffuse than 
for the two states of ZnO. Thus the bonding is much as expected and there are no 
changes in the bonding between the three states considered in this work. A CSOV 
analysis at the self-consistent-field level shows only small changes in the importance 
of Zn to 0 donation, 0 to Zn donation and of the strength of the cr bond. 

In the 3 n state, there is single cr bond and the bonding is enhanced by 0 2p7r 
donation to Zn. In ZnO - , the a bonding still occurs but with an additional electron 
in the oxygen 2 p 7 r orbital, the polarization of the 9 cr is much smaller, thus the cr 
bonding is more covalent and less ionic. The extra electron in the 4 tt orbital results 
in more 0 7 r donation to the Zn. Thus the bonding in ZnO - is stronger than in the 
3 n state, and the relative D 0 , r e , and u> e values are consistent with this. 

The bonding in the 1 E + state can be viewed as arising from Zn + 2 5(4s 1 ) + O - 
2 P(2pcr 1 2pTC 4 ) or Zn 1 S + 0 1 D(2pa°2pir 4 ). The population suggest that the bonding 
is derived from about equal mixtures of these two asymptotes. In the ionic limit, 
a bond forms between the Zn 4s and O 2 per orbitals and no promotion energy is 
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required to form hybrid orbitals. The 0 27 t donation is larger than for the II state 
since 0“ should donate more electrons than 0. In the 1 S+ l D limit, the Zn donates 
45 electron density to the empty 0 2 pa and the 0 2p7r orbital donates electrons to the 
Zn. The bonding is expected to be reasonably strong, and like the ionic asymptote 
no Zn hybridization is required. Thus the bonding in ZnO relative to these 

two asymptotes, is stronger than in ZnO relative to the Zn+0 asymptote. The 
D 0 of ZnO reported in this and previous work 3 ’ 5 is relative to Zn l S and 0 3 P, 
buth this ignores the promotion energy to reach the asymptote to which the 1 E + 
state formally dissociates. The 0 1 D state is 1.97 eV above 3 P (note the Zn + +0 
asymptote is even higher at 7.93 eV above the ground state asymptote). Therefore if 
the bonding in ZnO 1 S + is compared to its asymptote, the D 0 and oj e decrease and 
r e increases when ZnO" is formed by adding an electron to an antibonding orbital. 

IV. CONCLUSIONS 

The computed CCSD(T) spectroscopic constants are in good agreement with ex- 
periment. The calculations show that the 3 II state is not as low-lying as indicated by 
previous calculations. The ZnO EA is significantly larger than that of 0, despite the 
extra electron going into an antibonding orbital. This is consistent with ZnO 
formally dissociating to Zn 1 S and 0 l D. Thus relative to the appropriate asymptote, 
the Do and u> e of ZnO" are smaller and the r e longer than in the state of ZnO, 

which is consistent with the electron going into an antibonding orbital. Relative to 
the 3 n state, the additional electron increases the 0 to Zn 7r donation, and as a result, 
ZnO" is more strongly bound, has a larger co e and shorter bond length than the 3 II 
state of ZnO. The B3LYP results are in reasonable agreement with the CCSD(T), 
but important differences are observed. 
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TABLES 


TABLE I. Summary of computed results. 


Molecule 

»*e(A) 

w“(cm x ) 

w e Xe(cm x ) 

A>(eV) 

T e (eV) 

EA b (eV) 

CCSD(T) 







ZnO 1 E+ 

1.719 

727.2 

5.83 

1.63 


1.97(2.03) 

ZnO 3 n 

1.857 

566.6 

4.36 

1.38 

0.26 


ZnO” 2 S+ 

1.764 

664.4 

3.94 

2.20 



B3LYP 







ZnO X S+ 

1.713 

741 


1.21 


2.29(2.14) 

ZnO 3 n 

1.890 

509 


1.31 

-0.08 


ZnO- 2 S+ 

1.780 

592 


1.95 



CPF 5 







ZnO x £+ 

1.771 

646 


1.16 



ZnO 3 n 

1.866 

587 



0.03 



a For 64 Zn and 16 0. 

b The value in parentheses is corrected for the error in the 0 EA. 
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TABLE II. Summary of MCPF populations. 


State 9cr lOcr 4x 



Zn 

0 

Zn 

0 

Zn 

O 

Zn 

0 

ZnO 1 £ + ° 

0.93 

1.06 



0.44 

3.42 

+0.48 

-0.48 

ZnO 3 n 

0.55 

1.40 

0.84 

0.15 

0.10 

2.85 

+0.46 

-0.46 

ZnO” 2 E+ 

1.28 

0.63 

0.63 

0.37 

0.34 

3.54 

-0.27 

-0.73 


“ The X S + state has two important configurations and the populations from the 9(7 
and its correlating orbital are summed. 
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FIGURES 


FIG. 1. The MCPF 9 a natural orbital for the 3 II state of ZnO. 


FIG. 2. The MCPF lOcr natural orbital for the 3 II state of ZnO. 


FIG. 3. The MCPF natural orbital for the 3 II state of ZnO. 


FIG. 4. The MCPF 9a natural orbital for the state of ZnO. 


FIG. 5. The MCPF 47 t natural orbital for the X 1 S + state of ZnO. 


FIG. 6. The MCPF 9a natural orbital for the X 2 £ + state of ZnO . 


FIG. 7. The MCPF 10<r natural orbital for the X 2 £+ state of ZnO”. 


FIG. 8. The MCPF 47 r natural orbital for the X 2 £ + state of ZnO . 
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